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Abstract. To describe the interaction of molecular vibrations with electrons at a 
quantum dot contacted to metallic leads, we extend an analytical approach that 
we previously developed for the many-polaron problem. Our scheme is based on 
an incomplete variational Lang-Firsov transformation, combined with a perturbative 
calculation of the electron-phonon self-energy in the framework of generalised 
Matsubara functions. This allows us to describe the system at weak to strong coupling 
and intermediate to large phonon frequencies. We present results for the quantum 
dot spectral function and for the kinetic coefficient that characterises the electron 
transport through the dot. With these results we critically examine the strengths and 
limitations of our approach, and discuss the properties of the molecular quantum dot 
in the context of polaron physics. We place particular emphasis on the importance of 
corrections to the concept of an antiadiabatic dot polaron suggested by the complete 
Lang-Firsov transformation. 
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1. Introduction 



Recent advances in nanotechnology have stimulated great interest in the basic 
mechanisms of transport through molecular junctions (Chen et al. 1990, Park et al. 
2002, Reichert et al. 2002, Kubatkin et al. 2003, Park 2007, Cuniberti et al. 2005). In 
such devices the central element can be a single organic molecule or a suspended carbon 
nanotube, which may be thought of as a quantum dot contacted to metallic leads that act 
as macroscopic charge reservoirs. Transport through such a quantum dot is determined 
by energy level quantisation as well as electronic correlations and electron-phonon (EP) 
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interaction (Alexandrov & Bratkovsky 2003, Mitra et al. 2004, Takei et al. 2005, Galperin 
et al. 2007, Fehske et al. 2008). 

Vibrations of a molecular quantum dot are local excitations of substantial energy, 
which are represented by optical phonons. Their frequency is comparable to the transfer 
integral or kinetic energy of electrons (Nufiez Regueiro et al. 2007). Therefore, the 
mobility of electrons is significantly modified by the influence of molecular vibrations. In 
this respect, a molecular quantum dot resembles the situation in a crystalline structure, 
where the coupling between vibrations and electrons may lead to the formation of (small) 
polarons, as studied in the context of Holstein's molecular crystal model (Holstein 1959a, 
Holstein 1959&). A Holstein polaron is an electron dressed by a phonon cloud. Since the 
polaron must carry the accompanying deformation through the lattice, the mobility of 
Holstein polarons can be renormalised by several orders of magnitude in comparison to 
the free electronic excitation (Lang & Firsov 1962, Wellein & Fehske 1998). While the 
physics of Holstein polarons in a perfect crystal at low temperature and small density 
is by now well understood (see e.g. the review (Fehske & Trugman 2007)), there is less 
understanding if the periodicity of the crystal is altered, e.g. by impurities (Mishchenko 
et al. 2009, Alvermann & Fehske 2008) or disorder (Bronold & Fehske 2002, Bronold 
et al. 2004), in anisotropic materials (Alvermann et al. 2008, Emin 1986), or is absent 
for complicated geometries. 

For a molecular quantum dot, translational symmetry is broken from the outset, and 
EP coupling is relevant only in a small part of the entire system. The electron current 
through a deformable quantum dot was found to depend significantly on the local EP 
coupling (Flensberg 2003, Nunez Regueiro et al. 2007, Zazunov & Martin 2007, Mitra 
et al. 2004, Takei et al. 2005, Hohenadler & Fehske 2007). In order to understand the 
basic transport mechanisms in such devices, appropriate theoretical models have to be 
studied. The most simple model corresponds to a modified Fano- Anderson model, where 
a vibrating quantum dot replaces the static impurity. Then the current is determined 
by the dot spectral function (Meir & Wingreen 1992). The spectral function accounts 
for the leads, as well as for the influence of EP coupling. In particular, it determines 
the charge carrier population of the dot, and its value close to the Fermi energy of 
the leads determines the number of electrons contributing to the current. Since, with 
increasing EP interaction, spectral weight is transferred to lower energies, the charge 
carrier population of the dot increases. For the current, on the other hand, a reduction 
is expected since the spectral weight at the Fermi energy decreases. 

In the present manuscript we will address electron transport through a deformable 
molecule within an approximate description, which we previously developed for Holstein 
polarons at finite density (Loos, Hohenadler & Fehske 2006, Loos, Hohenadler, 
Alvermann & Fehske 2006, Loos et al. 2007). It accounts for renormalisation of transport 
and inelastic processes, and Pauli blocking. Higher order many-particle processes, 
namely the further excitation of electron-hole pairs and subsequent evolution of many- 
particle correlations, are not included. The current presentation, therefore, should 
be considered as an important but intermediate step towards a complete description. 
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Figure 1. (Colour online) Sketch of a molecular quantum dot with vibrational 
frequency wq between two metallic leads. 



A particular feature of our approach is that it interpolates between weak and strong 
coupling using an incomplete variational Lang-Firsov transformation. As a consequence 
it describes polaronic effects without being restricted to the antiadiabatic strong 
coupling regime. We introduce our approach here for the current in linear response, 
where the kinetic coefficient is obtained from the dot spectral function at equilibrium. 
Subsequent work will address the current at finite voltage bias. 

The paper is organised as follows. In section 2.1 we introduce the model 
Hamiltonian, and describe the variational Lang-Firsov transformation. In section 2.2 
we derive the expressions and the iterative calculation scheme for the self-energy, which 
depends on the variational parameter of the incomplete Lang-Firsov transformation. 
This parameter is obtained from minimisation of the energy, which we express as 
the expectation value of the Hamiltonian within the approximation used. From the 
spectral function, the kinetic coefficient is obtained in section 2.4. Section 3 discusses 
the numerical results, and we conclude in section 4. 

2. Theoretical Approach 



The paradigmatic example of a vibrating quantum dot is provided by a molecule 
sandwiched between two metallic leads (see figure 1). Such a system can be described 
by the Hamiltonian 



Here, the Eka (for k = 1, . . . , A^) give the energies of non- interacting electrons in the 
left and right lead a = l,r, and (Cf^^) are the corresponding creation (destruction) 
operators of free fermions in the A^ lead states. The leads will later be specified by their 
density of states q{^) = j^J^k^l^ ~ ^k]] the population of the leads is determined by 
the chemical potentials The quantum dot is represented by a single energy level A, 
with Fermi operators d^'^^ and the chemical potential /!=(/!; + fir)/'^- The term oc td 



2.1. Model 




(1) 
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allows for lead-dot particle transfer; a possible /c-dependence of the dot-lead-coupling 
can be absorbed in q{C)- An electron at the quantum dot interacts via a Holstein-type 
coupling with a local (intra-molecular) vibrational mode; g denotes the dimensionless 
EP coupling constant, and ujq the frequency of the optical phonons created (annihilated) 
by 6t ih). 

The quantum dot responds to the presence of an electron with a finite deformation. 
For sufficiently large phonon frequency cuq, the strength of the deformation depends 
only on the momentary occupancy of the dot. This is in analogy to Holstein's small 
polaron theory, where a lattice deformation in the vicinity of the electron accompanies 
the electron motion. To describe this effect, we apply a generalised Lang-Firsov 
displacement transformation (Lang & Firsov 1962) with parameter 7 G [0, 1], 

U ^^g{h^-b)dU ^ for g = ig. (2) 

After this transformation, the original electron and phonon operators are given as 

d = e9{b^-^)d, b = b + gdU. (3) 
The transformed Hamiltonian H = WHU reads 

k,a k,a 

+ {A-fi)dU -CddU + uJob'^b , (4) 



where 



and 



C, = -=t,e-^^''-'^ , C, = gujoa " l){b^ + b) , (5) 
V iV 



A = A - ep7(2 - 7) , with Ep = g'^uo . (6) 

As the parameter 7 of the Lang-Firsov transformation grows from 7 = to 7 = 1 
it accounts for the transition between the weak-coupling and strong-coupling regimes. 
The value of 7 will be later determined from minimisation of the energy. Only for very 
strong coupling and large phonon frequency, the value 7 = 1 is approached. Then, 
the canonical transformation (2) eliminates the direct coupling term between the new 
fermion and shifted boson operators at the price of introducing a boson modified transfer 
term between quantum dot and leads. This corresponds to the strong-coupling limit 
of polaron theory, where the new Fermi operators d would represent small polarons 
in the deformable lattice. For our problem, 'polaron formation' at the quantum dot 
mainly results in lowering of the dot energy level by the polaron shift Sp, and in an 
exponential reduction of the effective dot-lead transfer id = t^e"^^/^. Note that the 
variation of 7 throughout the parameter regime is important to describe the system 
away from the strong-coupling limit. The use of the 7-dependent variational Lang- 
Firsov transformation is an essential feature of our description. 
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2.2. Single-particle properties: quantum- dot spectral function 

We first determine the retarded Green function G^^ of the quantum dot, which is 
represented by the operators d^^'^ in the transformed Hamiltonian (4). The Green 
function is calculated within perturbation theory up to second order in the interaction 
coefficients (5), starting from the Lang-Firsov transformed Hamiltonian (4). Since the 
parameter 7 is assigned variationally, this treatment exceeds standard weak-coupling 
or strong-coupling perturbation theory, which starts either from the untransformed 
Hamiltonian (corresponding to 7 = 0) or the fully transformed Hamiltonian (7 = 1). 
The combination of perturbation theory with an incomplete variational Lang-Firsov 
transformation provides meaningful results also away from these limiting cases. 

Our calculation is based on the equations of motion for the generalised 
temperature Green functions (Kadanoff & Baym 1962), adapted to systems with EP 
interaction (Bruevich & Tyablikov 1962, Schnakenberg 1966). Accordingly we define 

G,,(ri,r2;{y}) = -^(T.d(ri)dt(r2)5), (7) 

and in an analogous way Gcd-ka, GdcM^ and Gcc;ka- The mean value and (imaginary) time 
dependencies in (7) are determined by H with fi = fii = fi^, the equilibrium chemical 
potential of the system. Moreover we set 

5 = exp I - dr V, (r)Q (r) + V, (r)Cl(r) + V, {t)C,{t) | , (8) 

where (5 is the inverse temperature, and the classical variables , , and are 
introduced ClS db purely formal device. 

We set up the equations of motion for the Green functions using the following 
matrix notation 

r dr'Gi(ri, t'; {V})G2{t' , r,- {V]) ^ G,{n, r'; {V])oG2{r\ r,:, {V}) .(9) 
Jo 

If Gi, G2 satisfy the relation 

Gi(ri, t'- {V}) o G2(r', r^, {V}) = 6[n - r^] , (10) 

they are called inverse functions of each other. In particular, the inverse functions to 
the zeroth-order Green functions G^^}{ti,T2) and Gf}f,^{Ti,T2) are given as 

d 



and 



Gir(n,r2) 



(A-/i) 



5[ri-r2] (11) 



d 

{Eka - /i) 



r^, 5[ri-r2], (12) 

respectively. By functional derivation with respect to the auxiliary fields {V} we find a 
set of coupled equations, 

Gf,-\n, r) o G,,{t\ r,- {V]) = 6[n - r,] 
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- ^.(ri, {V})GM{ri, ra; {V}) + — ^G,,(ri, ra; {V}) 
-Y.^t{ri,{V})G,,,Uri,r2;{V}) 
5 

mfn] 



k,a 



+ EAT77:;rxG'c.;'^«(^i.^2;m), (13) 



Gftkl(r,y)oG^,,Ur',r,-{V}) = 

- CKn, {V})G,,in, T,; {V}) + ^_G,,(ri, r,; {V}) , (14) 

With C,{t,{V}) = ^{%G,{t)S), Ct{T,{V}) = ^{%Ct{T)S), and Cj{r,{V}) = 
j^{%GKr)S). 

In order to solve this system of equations, we multiply (14) by G^J./^^ from the 
left and substitute the resulting expression for Gcd-M in (13). Then equation (13) is 
multiplied by G~^l from the right. The resulting equation for Gdd is converted to an 
equation for the self-energy Sd^, introduced by 

Gtd{n,T,-{V})=G'';;}-\n,T,)-T.dd{r,,r,-{V}). (15) 
By use of the functional differentiation rules 

5G o G-^ = -G o 5G-^ = G o 5G = Go5T.oG, (16) 
the self-energy S becomes 

S,,(ri, ra; {V}) = -Cd{rv, {^})5[ri - rj (17) 

+ 5^a(n; {V})G'^X{n,r,)Cl{r,- {V}) 

k,a 

- Z^'^cc,ka\^l^^2) T777-\ +Gdd[ri,r , {V\) o , . 



- J: Qin- {V})GZairu r") o Gdd{r\ r'- {V}) o ^^^^(^^^-^^ 

- E ^S..(n, r")a(-' ; {V}) o G..(r", r'; {r}) o ^""''^jll^^'f^ 

-|- terms with products of functional derivatives of T^dd 

+ terms with second functional derivatives of S^^. 

Within our iterative scheme, the terms on the r.h.s. of equation (17) without functional 
derivatives of Udd are taken in the first step as S^^''(ri, T2] {V})- Explicitly, we have 

S2(n,r2; {V}) = -Cdin; {V})6[n - r,] (18) 
+ E A (^i; {V})G%airu r2)Cl{r,; {V}) 

k,a 

'^'^ -'^^^{rrCAn)Cl{T,)S)-CAn-,{V})Cl{r,-,{V}) 



+ E^cta(n,r2) 



k.a 
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In the second step, to obtain the first approximation t!'^} has to be inserted into 
the functional derivatives of T^dd on the r.h.s. of (17). But, if we confine ourselves to 
terms up to second order in the interaction coefficients, only the first four terms on the 
r.h.s. are relevant. In this approximation, the self-energy Edd(ri,r2) = S(id(ri,r2; {0}), 
determining the temperature Green function Gdd(ri,r2) = Gdd{Ti,'T2', {0}), is given as 

S,,(ri,r2) = S2(ri,r2; {0}) + G,,in,T2){%Cd{n)Cd{T2)) , (19) 

where ^.(r; {0}) = and C,{r-{0}) = [C/(r; {0})]* = -^t,exp{-1^2 ^oth(i/3a;o)}. 
The correlation functions of the interaction coefficients occurring in (19) have been 
calculated previously for the generalised Lang-Firsov transformation in the Holstein 
model (Fehske et al. 1997). Converting (19) to the equation for the Fourier transform 
of the self-energy and expressing the Fourier transform of Gdd, Gddi^^u), by means 
of the spectral function Add, the summation over the bosonic Matsubara frequencies 
ujn = 2n7r / (3 (being the difference of two fermionic Matsubara frequencies) can be carried 
out, and we obtain, within low-temperature approximation Pujq ^ 1, 

i:^P^.)4^1e-^-E ; ,20) 

k,a 

V^g_g2 (^^^2 f np^Eka - n) _^ 1 - nF{,Eka - /i) 

^ s\ N ^ \iu;^ - {Eka - fJ') + suo iu^ - {Eka - 1^) - slvq 



/ + 00 / 
du'Addioo') - 



np{uj') ^ 1 — ripioj') 



^ lUJy — U + UJq lUJi, — UJ — UJq^ 

with the fermionic Matsubara frequencies uj^ = (2z/ + l)7r//3 and the Fermi function 
rip = (e^'^ + 1)"^. Analytical continuation icu^ ^ = + i5 in the upper complex 
half-plane then gives the retarded Green function 

Gf.(^) = - \ y (21) 

and the related spectral function 

Add{oo) = --lmGf,{uj + iO+). (22) 

71 

The r.h.s. of equation (22) is determined by the real and imaginary parts of Erfrf(ix' + iO"'"), 
which we obtain from (20) employing = V\ — 'm5{x). Moreover, we transform 

the /^-summation into an integration over the band energy ^ of the leads, ^ G [—W, W], 
using the lead density of states ^(0- Assuming further the right and the left leads to be 
identical the summation over a gives simply a factor 2. In the end, we work in the limit 
T ^ 0, when the Fermi function becomes the Heaviside function, nF{'jj) = 9(— tu), and 
obtain 

lmJ:ddiuj) = -27itle-3'g{uj + fx) / d^ - (^ - /i)] (23) 

J-w 

-27rt2e^^' V^^i^(cu + /i + swo) / 6[uj - - fi) + sluq] 
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- TT [(1 - 7)5-^^0]^ {Add{'^ + ujq)<c){-uj - ujo) + AddiuJ - ioq)<c){uj - LOq)} 



2tle--^' 



S>1 



W 1 

de^^(0 — ^ — r 

■w ^-[i-N 

' def^(0 ^^-TT 

_w uj-[i- 11) + sujQ 

+ / ^iQiO ^ ^ 



(24) 



duj'- 



g UJ — UJq — UJ 



Note that in accordance with our second order approach, the spectral functions Addif^) 
occurring in (23), (24) have to be determined by equations (21), (22), using t!{^{uj + \Q^) 
for the self-energy. 

For = 0, we are faced with the well-known problem of electron localisation at an 
impurity. Then, if a solution uj of 

u; = A + Re ^^^(tD - /i) (25) 

exists outside the interval [—W, W], the spectral function Addiyj) exhibits a single-peak 
structure, 

Addif^) = z{u)5[uj -{u - /i)] , (26) 
reflecting electron localisation at the quantum dot. Here, 

1 



ZiUJ 



l + 2tlV 



w 



w 



(27) 



According to the first term on the r.h.s. of (23), we find \m.'Ldd{'^) 7^ for G 
[—W — /i, — /i], leading to an incoherent single-particle spectrum in this interval. 

For non-zero coupling between the electron and the local vibrational mode at the 
dot, we have ImSrfd(t^) 7^ everywhere provided that ojq < W — Then, if (7 7^ 0, 
the spectral function Addif^) exhibits no coherent contribution, and is given, for all uj 
and 7 G [0, 1], by a purely incoherent spectrum 

1 ImSdd(u;) 



Addit^) 



u -{A- n)-Re T,ddiuj) + [Im ^^^(u;)] 



(2^ 



-(S+I)(0o 



-SWo 



-w 



w 



SWo 



(8+1)0)0 



W+ 1^ W+ |x 



W-|i 



-► CO 



W-|x 



Figure 2. Shaded regions indicate the intervals of non-zero contributions to Im'Eddiijj) 
from the processes determined by the first two terms on the r.h.s. of equation (23). 
Note that all intervals overlap if cuq < W — Shown is the case uiq > W, /x < 0. 
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On the other hand, if the latter condition is not fulfilled because —ujq < —W — fi 
OT W — fi < Uq, then the contributions to ImT^ddi^^) from the phonon processes 
vanish in certain cj-intervals (see figure 2). In particular, this happens for c^; < 0, 
uj e [— (s + l)a;o, —suq — W — fi] and for > 0, a; G [W — /i + sujq, (s + l)u;o] with 
s > 0, respectively. Poles in the Green function occur if equation (25) has a real- valued 
solution in these intervals. A small width of these peaks can arise from the third term 
on the r.h.s. of equation (23) if 7 differs appreciable from unity. 



2.3. Determination of 'y 

In order to fix the variational parameter 7 self-consistently, we minimise the ground- 
state expectation value E = (H) with respect to 7. We factorise the statistical 
averages with respect to phonon and polaron variables, i.e. (Cdd^d) ~ {Cd){d'^d) and 
{Ctd^Cj^a) ~ {Ct){d^Cj^a) and obtain 

E = J2{E,^ - n){clcj - J2 {{Ct){d'cJ + (C/)(cU) 

k,a k,a 

+ (A - fi){dU) - (Cd) (dU) + uoib^b) . (29) 

For the transformed Hamiltonian, (b^b) as well as (Cd) is zero for T ^ in our second 
order approach. The remaining expectation values can be expressed as 



(dU) = / duj'Addiuj')nFiuj'), (30) 
J —00 

{d^Cka) = Gcd;kaiTi,T2) _= / dJ Upiuj') ImG ^d-^ ka{uj') , (31) 



oo 
oo 



1 1"^ 

{cla^ka) = Gcc;ka{Tl,T2) _ = / duj' np^u') Im G cc; kai^j') . (32) 

Based on the equations of motion (equation (14) for Gcd-ka and an analogous one for 
Gcc;ka), the Grceu functions are determined to first order as 

G.dM^) ^ - G%S^){Cl)Gddiz) , (33) 

Gcc;.a(.) ^ G^Xi') + {Gl) {G, ) [G^X(')] 'Gd,{z) , (34) 

which is consistent with the factorisation of bosonic and polaronic averages in (29). Via 
Gdd, the expressions (33) and (34) depend on the spectral function Add- We identify the 
variational parameter 7 from the minimum of -^(7). 



2.4. Two-particle properties: kinetic coefficient 

In this section we calculate the current of the system (1) caused by a small potential 
difference between the leads. To this end we deduce the expression for the kinetic 
coefficient in terms of the quantum-dot spectral function, consistently with the 
approximations assumed in the derivation of Add{i^)- In accordance with general linear 
response theory (Zubarev 1971), the external perturbation Ht = —aF{t) coupled to the 
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system variable a induces the change (d) of (a), whose Fourier transform (d)^ is related 
to the Fourier component F{uj) via the kinetic coefficient L^uj) as 

(d), = L{uj)F{uj) . (35) 

In particular, linear response theory gives for L = L{uj 0) the expression 



lim 



-— Im ((d|d)), 

UJ 



(36) 



where the symbol {{■\-))uj denotes the retarded (commutator) Green function in 
frequency representation. In our case, the quantity of interest is the current between 
the lead a and the dot. By the continuity equation, the operator for this current is 

Ja = -eNa , (37) 
where iV„ = J2k Aa^ka and 

Na = i[H, Na] = -i-^ td J2 {d^^ka ' clJ^ (38) 



with the electron d operators defined in (3). With a = iVa in the above general linear 
response formula, we express the kinetic coefficient, obtained from variation of the 
chemical potential /Xq = /i + A/i^ in the first term of the Hamiltonian, through the 
retarded Green function {{Na\Na)). 

To obtain this Green function, let us consider the corresponding Matsubara Green 
function 

GAn, r[) = - (T,iV„(ri)iV,(rO) (39) 
which is related to the two-particle Green function by 

Gj{n,T[) = G{r2,r^,n,T[) _ , (40) 



where 

<5(r2,r2;ri,r{) = -^^1^ \(TrCka{T2)ck'a{T2)d\T[)d\Ti] 

k,k' 

- (TrCka{r2)d{T2)cl,^{T[)d\Ti] 

- (%d{T2)Ck'a{r2)d\T[)cl^{Ti) ^ 

^T.d(r2)d(r^)4,,(r0cL(ri))j (41) 

The mean values of all the time-ordered products on the r.h.s. of (41) are basically two- 
particle Green functions G{2, 2'; 1, 1') (Rickayzen 1981), which may be approximated by 
one-particle Green functions according to 

G(2, 2'; 1,1')^ G(2, 1)G(2', 1') - G(2, lOG(2', 1) , (42) 

if vertex corrections due to phonon-mediated electron-electron scattering are neglected. 
In this way, the terms on the r.h.s. of (41) turn out to be proportional to products of 
Gcd^ ^dc' Gcc: and G^j. Substituting d from (3), averaging over the oscillator variables 
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and inserting the expression (33) for Gcd and G^^ = G*^, it becomes obvious that the 
terms containing the latter "mixed" Green functions are of higher order in |(Ct)|. Hence 
we get as the leading order result 



k 



(43) 



Here we have used the unperturbed Green function G'^^J (see equation (12)) for the 
electrons in the leads and introduced the notation Gdd = G^^. 

Inserting (43) into (40) and performing a Fourier transformation, the latter equation 
becomes 



GjiiiVn) = J^tdYlj^Yl + iWn)Grfd(iw^) 



(44) 



with bosonic Matsubara frequencies Un- We now express Gdd by the electronic spectral 
function Add{^^) 



Gddi^^p 



make use of 



1 



(45) 
(46) 



perform the Matsubara summation over the fermionic frequencies ujy, and obtain 



1 f p + oo 

Gj{iu;^) = -tlJ2\ / do;' A 

7 K.J —CO 



k 

+ CXJ 



i.upiuj ) - npiEka - At) 
ddl^ )-■ : — ; — 7^ r 



du'AddiuJ , . 

^ lUn -UJ' + [Eka - /i) 



(47) 



The analytical continuation of (47), iu^ ^ u = 00 + 16, gives the retarded Green function 

{{Na\Na))^ = Gj{uj + iO+), (48) 



leading to 



Im((iV,|iV,,)), 



N 



tl ^ ^^Add{Eka - fi-uj) [npiEka - fi-uj) - npiEka - /U)] 



+Add{Eka - fi + uj) [ripiEka - At) - np^Eka - fi + uj)] 



(49) 



Assuming identical leads a = l,r, this result is of course independent of a. Finally, 
according to the definition (36), we have to perform the limit — > 0: 



L = lim 

a;->0 



Im((iVjiV,)), 

00 

w 



211 
iV 



tl J2 ^dd{Ek - /i) {-n'p{Ek - At)) 



2vrt^ / de giOAddi^ - fi) {-n'A^ - ft)) 
-w 



(50) 
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Then, for T — » 0, (— — jf)) = — fi], so that 
L = 27itlg{fi)AM{0). 



(51) 



The general relation between the electronic spectral function Add{uj) needed here and 
the polaronic spectral function A^iuj) determined in the preceding section, was derived 
in Ref. (Loos, Hohenadler & Fehske 2006). Using equation (40) of this work, for u = 
and T ^ 0, it simply follows that 



Consequently, the kinetic coefficient is determined by the value of the polaronic spectral 
function Add{^^) at the lead Fermi level, multiplied with the renormalisation factor e~^^. 
In particular, changing g ot A leads to a shift of A^di^^) with respect to the lead Fermi 
level and therefore changes the value of Add{0), as will be shown by the numerical 
calculation in the following section. 

3. Numerical results 

While the derivation of the equations in section 2 is completely general, we consider 
in the following the case of a single quantum dot between semi- infinite ID leads, as 
sketched in figure 1. Accordingly, the density of states of the leads is given by 



with the half-bandwidth W = 2t. We fix t = 1 from here on. 

The numerical computation of the Green function Gddi^^), equation (21), and the 
corresponding spectral function Add{^^), equations (22), (28), is performed by evaluation 
of equation (20) for energies u = u + iS slightly above the real axis. A small choice of 
S > avoids problems arising from the simultaneous treatment of poles and incoherent 
parts in the Green function. For our computations, we used 6 ^ 10~^. Alternatively, one 
might directly evaluate equations (23), (24) which are given for real u, but numerical 
inaccuracies in the calculation of the principal value integrals tend to degrade the 
computation. From the spectral function, we obtain the kinetic coefficient L using 
equations (51), (52). 

It is noteworthy that, restricting ourselves to second order perturbation, the 
approximations to the Green functions preserve important sum-rules for the spectral 
functions, especially for the integrated weight (Koch 2009). This remains true for finite 
5 > in the numerical computation (while the numerically unfavourable evaluation 
of equations (23), (24) led to the artificial drop of total spectral weight reported in 
Ref. (Loos, Hohenadler & Fehske 2006)). 

In the discussion of the numerical results we start with important limiting cases. 



AddiO) = e-'s'AddiO) ■ 



(52) 




(53) 
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Figure 3. (Colour online) All results for = and /i = 0. Left panel: Dot spectral 
function Add{'^) for A = and several td- Right panel: Kinetic coefficient L as a 
function of A for several = 1. 



3.1. Non-interacting case 

For vanishing EP coupling [Sp = 0) the problem reduces to that of an impurity in a 
ID chain. The spectral function is then obtained exactly by our calculation (left panel 
in figure 3). For t^ = 1, the spectral function of the translationally invariant ID chain 
is recovered. For smaller t^, the spectral function develops a pronounced maximum at 
uj = 0, which evolves into a 6-peak in the limit — 0. 

For the kinetic coefficient, we observe in figure 3 (right panel) two effects that will 
be important later in our discussion of the interacting case. First, for dot energy A 7^ 
scattering off the dot impurity leads to reduction of L compared to the case A = with 
minimal scattering. Consequently, L is maximal for A = and shrinks monotonically 
with growing |A|. Second, a reduction of leads to a reduction of L. Moreover the 
variation with A becomes more pronounced as electrons become more susceptible to 
scattering off the dot. In equations (23), (24) we see that for smaller hybridisation 
td^iO ^he broadening of dot levels due to coupling to the continuum of lead states is 
reduced. In the limit t^ — > 0, the function L = L{A) becomes a 5-function at A = with 
weight oc t^. Note that L is independent of td for A = 0, which is however a peculiarity 
of the non-interacting case without damping of states close to the Fermi energy. 

3.2. Small phonon frequency 

In the following, we first discuss the results of our approach in the limits of small and 
large (next subsection) phonon frequencies. For the moment, we fix = t = 1 and 
/i = which corresponds to the half-filled band case for a translational invariant system 
with Ep = and A = 0. 

For small phonon frequency Uo/t = 0.1 (adiabatic regime), when the phononic 
timescale is much slower than the electronic timescale, we expect significant deviations 
from the behaviour described by the standard Lang-Firsov approach consisting of a 
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complete Lang-Firsov transformation and a subsequent average over the transformed 
phonon vacuum. Our approach is able to account for these deviations by the variational 
parameter 7. The deviation of 7 from unity is some measure of both adiabatic and 
weak-coupling corrections. 



1.5 




Figure 4. (Colour online) All results for A = 3, i = = l, = 0, and loq = 0.1. 
Upper panels: Total energy E (a) as a function of 7. for different £p (inset: Minimum 
of as a function of £p) and optimal parameter 7™^, kinetic coefficient L and particle 
density n at the dot (b) as a function of ep. Lower panels: Dot spectral function Addi'^) 
and integrated spectral weight S{ijo) for different £p (c)-(f). The insets in panels (c) 
and (f) show the imaginary part of the second-order self-energy Y,dd{^)- 

3.2.1. Repulsive dot In figure 4 we show results for the case uiq = 0.1, /i = 0, and a 
repulsive dot A = 3. Since the dot is repulsive, the particle density at the dot is small, 
n ^ 0.1, for small EP coupling Ep [see panel (b)]. At a critical coupling ~ 3.325, EP 
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interaction at the dot overcomes the repulsive potential, and a transition takes place to 
a situation with large n. This transition is accompanied by a jump of the variational 
parameter 7mm from a small value (< 0.02) to 1. This jump can be traced back to 
the behaviour of the total energy E a.s a. function of 7: If Sp increases, -^(7) develops 
two local minima [see panel (a)]. At Sp = Sp, the minimum at 7 = 1 becomes the new 
global minimum. Evidently, the sudden change of 7 reflects the formation of a strongly 
localised polaron at the quantum dot. Thereby the lead-dot transfer is almost completely 
suppressed and, in accordance with this picture, the kinetic coefficient (oc exp{ — (7^}) 
drops to zero at the transition point. This suppression of transport is well described 
by the complete Lang-Firsov transformation (although this basically is a non-adiabatic 
approach), mainly because we enter an extreme strong-coupling situation [g'^ = 5.77). 
Note that the observation of an extremely sharp polaron transition in the adiabatic 
regime for repulsive quantum dots is in accordance with recent exact diagonalisation 
results (Fehske et al. 2008, Alvermann & Fehske 2008). 

We next analyse the dot spectral function Add [see figure 4 (c)-(f)]. For Sp < Ep, 
we have g < 'jming'^ < 0.1154 and, calculating H^JJ within our second order scheme, 
the first term in equation (20) is basically proportional to the semi-elliptical density 
of states of the leads, while the second term ("phononic contribution") is insignificant. 
Then the resulting spectral function A^JJ, which has to be put into the third term of 
(20), describes a continuum of states, roughly in between -2 and 2, and a localised dot 
uo ~ 3.5 (in accordance with the result obtained for the Fano- Anderson model). 
Since the prefactor of the third term, (1 — '~imin)g^o ^ 0.57 for g < g'^, is rather large 
this term gives a significant contribution to the second order self-energy E^^. Thereby 
the localised peak in A^JJ becomes evident in S^^ [see inset of figure 4 (c)]. As a result 
the second order spectral function. Add, exhibits two sharp peak structures (localised 
states) above the continuum of states around u = 0. If the EP coupling increases these 
peaks become more and more separated. In order to analyse the spectral weight of the 
different signatures in Add, we have calculated the integrated spectral function 

Siuj) = I duj'Add{uj') . (54) 
^—00 

Figures 4 (c) and (d) show that for Sp < Sp the spectral weight mainly rests in the 
localised peak structures above the wide band. Hence the spectral weight of the 
current-carrying states at the Fermi energy /i = is reduced, and the kinetic coefficient 
L oc Addif^) is substantially lowered compared to the case A = 0. At = e^, 7mm 
jumps to 1, and the strong renormalisation arising from the complete Lang-Firsov- 
transformation results in a pronounced peak at negative energy at about A — e^, which 
now, however, is the signature of a quasi-localised polaronic dot state. The polaronic 
quasiparticle peak is accompanied by two side-bands (roughly of width 2W) shifted by 
i^p, which arise from the Poissonian distribution of phonons at the dot, with maximum 
at g"^ = Ep/ ujq phonons. States in these band are strongly damped due to the significant 
phononic admixture, as is evident in the imaginary part of the self-energy [see inset 
figure 4 (f)]. 
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3.2.2. The case A = For A < the quantum dot is attractive. For A = and = 
we have of course a translational invariant ID system, where /i = corresponds to the 
half-filled band case, i.e., n = 0.5. Such a "neutral" quantum dot becomes attractive 
for arbitrarily weak EP interaction. This is because the "effective" dot level is given by 

A A £^p0'mm(2 'ymin) • 

Consequently, in figure 5 (b) the particle density at the dot is larger than 0.5 for all 
Bp, and the dot spectral function has no pole at positive energies. At small EP coupling 
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Figure 5. (Colour online) All results for A = 0, t = = 1, /i = 0, and loq = 0.1. 
Upper panels: Total energy E (a) as a function of 7, for different Ep (inset: Minimum 
of £^ as a function of Ep) and optimal parameter ^mim kinetic coefficient L and particle 
density n at the dot (b) as a function of £p. Lower panels: Dot spectral function 
Addii^) and integrated spectral weight S{lo) for different Ep (c)-(f). The inset in panel 
(f) gives the imaginary part of the dot self-energy Y,dd{^)- 
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the spectral function is similar to that of of a ID tight-binding model. The weak EP 
interaction causes the spiky signatures separated by ujq from the upper and lower band 
edges [see figure 5 (c)]. In contrast to a repulsive dot, for which our methods correctly 
describes the transition from unbound to localised polaronic dot states, a sharp polaron 
transition cannot occur for a dot level A < 0. Nevertheless, we observe in figure 5 a 
transition signalled by the jump of to 1 with corresponding increase of n, decrease 
of L, and formation of a pronounced peak in ^^^(u;) at negative energies. The reason is 
again the change of the global minimum of £'(7), which has two local minima for larger 
Bp. Since for A = the interaction need not overcome a repulsive dot potential, the 
transition takes places at smaller ~ 1.125. Therefore, and in contrast to the previous 
case, no isolated quasiparticle peak in Addijjj) emerges at the transition, and the change 
of the spectral function is less dramatic. Above the transition the qualitative behaviour 
of the imaginary part of the self-energy [see inset of 5 (f)] is the same as for the repulsive 
quantum dot [cf. inset of 4 (f)]. Since Sp is smaller now, the maxima of the phonon 
contributions to S^d are less separated than in figure 4(f). 

In our approach the transition results from a jump in 7mm- As before, this might 
indicate the formation of a localised polaronic dot state. But we know from the various 
variational approaches to the polaron problem that such jumps often arise as artefacts of 
the variational ansatz (Fehske et al. 1994). For the Holstein polaron with EP interaction 
at each lattice site, no phase transition exists (Gerlach & Lowen 1991). Instead, a 
crossover between an almost free electron and a heavy polaron takes place. The crossover 
can however be very rapid for small phonon frequency (Alvermann et al. 2008). But 
we also know that, in contrast to the Holstein polaron problem, for a single electron 
at a vibrating quantum dot a true phase transition, from n = to finite n > 0, takes 
place (Mishchenko et al. 2009, Alvermann & Fehske 2008, Fehske et al. 2008). This 
phase transition becomes more pronounced for small ujq. The behaviour found here 
therefore does not contradict the essential physical mechanism in our situation. In 
principle, our approach mimics the sharp adiabatic polaron transition by the change of 
the parameter 7 of the (non-adiabatic) Lang-Firsov transformation. While the precise 
nature of the transition is only poorly described by this approximation, we still believe 
that the transition - or rapid crossover - itself is characteristic for the quantum dot at 
small uiQ. 

3. 3. Large phonon frequency 

For large phonon frequency ujq = 10, in the antiadiabatic regime, phonons adjust 
instantaneously to the electrons. Now our non-adiabatic variational Lang-Firsov 
approach perfectly matches the situation. We will see that the transitions found in the 
previous (adiabatic) cases will be replaced by smooth changes of the physical observables. 
We then note that the results obtained can be understood easily starting from the case 
without EP interaction. 
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3.3.1. Repulsive dot The fact that for large phonon frequency no transitions occur is 
most clearly seen for a repulsive barrier in figure 6: All quantities depend smoothly on 
Bp. The total energy -£^(7) has a unique minimum for all e^, which is the reason why no 
transition occurs. Note that 7 grows from 0.75 to 1, increased. For large ujq 

and Bp, the Lang-Firsov transformation implements the correct physical mechanisms. 
Nevertheless, at weak EP coupling, the deviation 7 < 1 indicates the importance of 
corrections to the complete Lang-Firsov transformation. 

The spectral functions in figure 6 show that, although no transition occurs, we 




-20 -15 -10 



Figure 6. (Colour online) All results for A = 3, i = = 1, M = and = 10. Upper 
panels: Total energy E (a) as a function of 7, for different £p (inset: Minimum of E as 
a function of Sp) and optimal parameter ^mim kinetic coefficient L and particle density 
n at the dot (b) as a function of Ep. Lower panels: Dot spectral function A^ii{u;) and 
integrated spectral weight S{uj) for different Ep (c)-(f). The inset in panel (f) gives the 
imaginary part of the dot self-energy T.dd{i.^)- 
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start with a peak in A^di'^) at positive energies for small Sp (< A), to end up with a 
polaronic quasiparticle signature at very strong EP coupling. At = 2, the peak enters 
the band of lead states from above, leading to an asymmetric deformation of the semi- 
elliptic band [see panel (c)]. A second absorption feature is separated by the phonon 
frequency, but carries almost no spectral weight. With increasing EP interaction, the 
effective dot level is lowered until a "neutral" dot evolves at about A = A — = 0, 
as can be seen from the ID tight-binding-model-like absorption in panel (d) (note that 
we here are in the weak EP interaction regime since = 0.3, which is the relevant 
coupling parameter in the anti-adiabatic region, is small). At large EP coupling, the 
polaronic peak appears at A < and acquires a spectral weight of nearly unity (see 
panel (f) for Ep = 8). Due to the large phonon frequency ujq > W, the phononic 
sidebands do not overlap in this case (in contrast to figure 4), and the spectral function 
and self-energy show the typical multi-band structure known from the antiadiabatic 
Holstein polaron. Most importantly, we now find intervals where ImE^^^ = between 
the non-overlapping phonon (side) bands. If the polaronic peak is located within such 
an intermediate range the quasiparticle cannot decay by (multi-) phonon absorption 
or emission processes [see panel (f)]. This means the polaronic dot state acquires in 
principle an infinite lifetime (in the limit of very large couplings and phonon frequency). 
Naturally, as Ci;o ^ cxd, we recover the behaviour of the impurity model, where a true 
bound state occurs (Alvermann & Fehske 2008, Mishchenko et al. 2009). 

Thus, for large phonon frequency, or whenever 7 is close to unity, we can understand 
most properties starting from the non-interacting case, if we take the interaction into 
account by renormalisation of the appropriate physical parameters. Inspection of 
equations (51), (52) shows that one central effect of interaction on the kinetic coefficient 
is the renormalisation of to an effective dot-lead hopping t^e"^^/^. The second central 
effect is the change in the dot density of states, which is to a large extent caused by 
lowering of the effective dot energy (below the value A without interaction) due to 
deformation of the quantum dot in the presence of electrons. 

The simple picture is valid only in the limit 7 = 1, when the dot energy is 
effectively lowered by —Sp, such that A = A — in equation (6), and the dot-lead 
hopping is effectively reduced by e~^^/^, such that ^^^(0) = e~^^ Add{0) in equation (52). 
The kinetic coefficient L then has properties analogous to the non-interacting case, 
with the appropriately renormalised parameters. We discussed above (section 3.1) the 
consequences for L resulting from a change of or A. For the curve shown in figure 6, 
it turns out that it can be indeed reproduced from the expression for L in the non- 
interacting case, evaluated with an effective dot energy A — Sp and effective dot lead 
hopping td = t^e"^^/^ replacing A, td- In particular, L is maximal for Sp = A (cf. 
figure 3). Note that away from the limit of large phonon frequency, whenever 7 <^ 1, 
different behaviour is found. Also the shape of the dot spectral function, and especially 
the value of Add{0), is modified in addition to simple renormalisation. Of course, and 
similar as for the polaron problem, the retardation of the EP interaction manifests itself 
most prominently at small to intermediate phonon frequency. 
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3.3.2. The case A = The behaviour for the attractive dot is similar to the previous 
case (see figure 7). Here, of course, a pronounced peak in A^dioj) occurs at negative 
energies for all Sp. Once again, all features can be understood starting from the non- 
interacting case with appropriate renormalisation, as explained above. Since A < for 
all En > 0, the kinetic coefficient has no maximum as a function of €„. 
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Figure 7. (Colour online) All results for A = 0, = — ^^'^ = 10. Upper 

panels: Total energy E (a) as a function of 7, for different Sp (inset: Minimum of E as 
a function of £p) and optimal parameter ^mim kinetic coefficient L and particle density 
n at the dot (b) as a function of Ep. Lower panels: Dot spectral function A^ii{u;) and 
integrated spectral weight S{uj) for different Sp (c)-(f). The inset in panel (f) gives the 
imaginary part of the dot self-energy Y.dd{^)- 



The simple picture given above takes into account only the renormalisation of t^ 
and A. It is important to keep in mind that both effects lead to a reduction of L. As a 
consequence, the change of the kinetic coefficient is not simply given by an exponential 
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behaviour oc (compare L to the dashed curve in figure 7). In the present case, the 
couphng strength is small in terms of the average number of phonons g'^ = Ep/uQ, for 
which g^ < 1, but large in terms of the shift of the dot energy Sp, which is of the order 
of the bandwidth W. Here, the reduction of L is mainly caused by this large shift. 

The opposite situation can occur for small phonon frequency, when g"^ is large 
already for small Sp. Then, however, the renormalisation of is not adequately described 
by an exponential factor e"^^/^. In the limit t^o ^ of small phonon frequency, 
g^ = Ep/uJo oo for any Sp > 0. If the exponential dependence oc e~^^ persisted, 
that would imply zero current even for tiny Sp, which is unphysical. A calculation with 
fixed 7 = 1 therefore overestimates the reduction of L for intermediate-to-small phonon 
frequencies. We discussed in the previous subsection how, in our treatment, variation 
of the parameter 7 accounts partially for this deviation, leading to 7 -C 1 away from 
the antiadiabatic strong-coupling limit. 

3.4- Intermediate phonon frequency 

For intermediate phonon frequencies the qualitative behaviour depends crucially on the 
value of 7, even if no transition occurs. From our previous discussion we know that 
both a positive A or a small uq favour a rapid, or even discontinuous, transition. For 
cuo = 1, we show in figure 8 (upper row) how a smooth crossover evolves into a sudden 
transition with increasing A. In contrast to the case of small phonon frequency ujq = 0.1, 
the kinetic coefficient L is a smooth function of Ep for A = 0. A transition in L occurs 
only for larger A. Increasing the phonon frequency to ujq = 3 [lower row, panel (c)] then 
leads again to a smooth crossover even at A = 3. 

Changing the phonon frequency, we should ask to which extent the renormalisation 
scenario given for the antiadiabatic case remains applicable. For ujq = 3 [panel (c)] we 
observe that L differs from the value obtained, as in the previous subsection, from the 
non- interacting case for renormalised A, (in particular the maximum of L occurs for 
> A), but although A = 3, the two curves match rather well. The situation changes 
for uq = 1 [panel (d)], where strong deviations occur already for A = (note that 
the dashed curve for A = 2 even misses the increase of L at smaller Sp). Evidently, 
the simple renormalisation scenario fails, as we expected. We can achieve much better 
agreement if we perform the same calculation but incorporate the parameter 7 taken 
from the upper left panel in figure 8 (the dashed curves would correspond to fixed 7 = 1). 
Small deviations remain for A = 2, since the full calculation includes damping of states, 
indicated by a finite imaginary part of the self-energy, which is not captured by the 
change of 7. 

It is now evident that the essential feature of our calculation is the self-consistent 
determination of the parameter 7. Once we know its value, we may get a good 
approximation already with a modified renormalisation argument which was originally 
constructed for the antiadiabatic limit. If, in contrast, we fix 7 = 1 we will miss the 
physics away from the limit of large phonon frequencies. The restricted use of the Lang- 
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Figure 8. (Colour online) All results for t = <d = 1, A* = 0. Upper row: For wo = 1, 
optimal 7min (a) and kinetic coefficient L (b) for various A as indicated. Lower row, 
panel (c): For A = 3 and wq = 3, optimal parameter 7mm, kinetic coefficient L, and 
particle density n at the dot, as a function of Ep. The dashed green curve shows L 
calculated for the non- interacting case with renormalised parameters (see text) . Lower 
row, panel (d): For ojo = 1 and two different A, kinetic coefficient L as a function of £p. 
The dashed curves show L calculated for the non-interacting case with renormalised 
parameters, but fixed 7 = 1. The dot-dashed curves have been obtained taking the 
parameter 7 from the upper left panel (see text). 



Firsov transformation for intermediate-to-small phonon frequencies is well known in the 
Holstein polaron literature. It is important to realize that this restriction applies also 
to the situation of a vibrating quantum dot. 

3. 5. Variation of the chemical potential 

So far all results were given for chemical potential /i = 0. A change of the chemical 
potential affects the kinetic coefficient in two ways. First, since in equation (51) the lead 
density of states q{^) and the dot spectral function A^di'^) are evaluated at the chemical 
potential, a change of n results in a change of L. Second, phonon emission/ absorption 
is possible only if free states are accessible after an electron changed its energy by iscjQ. 
Otherwise, EP interaction is suppressed by Pauli blocking. Therefore, the shape of 
Addijjj) itself does depend on /i in a true many-particle calculation as performed here. 
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Figure 9. (Colour online) All results for A = 0, f = = 1 and cjo = 1- Comparison 
of Add{^), S{lu), and the self-energy 'Sddi^^) (insets) for /i = [left column, panels (a) 
and (c)] and /i = —1.9 [right column, panels (b) and (d)], for weak coupling Sp = 0.5 



[upper row, panels (a) and (b)] and strong coupling Sp 



3.1 [lower row, panels (c) and 



Significant changes occur whenever ±u!q crosses the band edges (at about ±W — /i at 
weak couphng). 
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Figure 10. (Colour online) All results for A = 0, t = td = ^ and luq = 1. Kinetic 
coefficient L [panel (a)] and dot density of states n [panel (b)] as a function of Sp, for 
varying chemical potential /.j. 

This effect is evident in the spectral function Addii-^) in figure 9. At weak couphng 
{Sp = 0.5) the shape of Add{^^) is similar for half-filling (/i = 0) and small particle 
density (/i = —1.9), but small differences at the lower band edge are a first indication 
of the different behaviour at stronger coupling. There, for Sp = 3.1, the spectrum for 
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yU = is completely incoherent, with finite Im Tj^di'^)- Around c<j = we observe a valley 
in ImSrfrf(u;) of width 2ujq, which results from Pauli blocking of states in the vicinity 
of the Fermi energy (cf. the discussion in Ref. (Loos, Hohenadler & Fehske 2006)). 
Note that ImSrfrf(u;) ^ here even at the Fermi energy, since the self-energy contains 
the contribution from dot-lead transfer. For /i = —1.9, states below the phonon 
emission threshold, located ujq above the lower band edge, cannot emit a phonon (phonon 
absorption is suppressed at zero temperature). Electrons in these states are undamped, 
with infinite lifetime corresponding to \m'Lcid{'^) = 0. 

The interpretation of the behaviour of the kinetic coefficient L (see figure 10) relies 
on these two mechanisms. First, if fi decreases, the change in the density of states 
should reduce the value of L [compare the curves for /i = (solid line) and fi = —1.9 
(dot-dashed line)]. Also, the dot density of states n decreases. We note that for the 
non- interacting {sp = 0) ID case the changes in ^(^) and Add{uj) cancel by chance, and 
L is independent of fi. However, at stronger coupling, the different infiuence of Pauli 
blocking reverses this behaviour, and L is larger for smaller /i. This explains why the 
curve for /x = —1.9 crosses the curve for /x = in figure 10. 

3.6. Small dot-lead hopping (tunnel contacts) 

We have so far discussed the importance of the phonon frequency only in the situation 
t = td- On physical grounds it is the ratio c^o/td, instead of ujo/t, which should distinguish 
the adiabatic from the antiadiabatic regime. 

In figure 11 (a) we show, for intermediate phonon frequency ujq = 1 and A = 0, the 
change of behaviour as td is reduced by one order of magnitude. The kinetic coefficient 
L decreases with t^ (see upper panels). Note that if we reduce td at fixed A = ep/{2td) 
in figure 11 (a), we de facto reduce the coupling e^, but nevertheless L decreases. This 
implies that the effect of smaller td on L dominates over the possible increase of L for 
smaller coupling. We compare L again with the value obtained in the renormalisation 
scenario, with a calculation as in the non-interacting case but with renormalised 
parameters A, t^ (see dashed lines). For t^ = t = 1, and c^o/td = 1, both curves disagree. 
We already discussed above that this is a consequence of adiabatic corrections, which 
become important at intermediate-to-small phonon frequency. However, if we reduce 
tfi to 0.1 while keeping ujq = 1 fixed, and thereby increase uj^/td to 10, both curves 
match. Apparently, we enter the antiadiabatic regime by sufficient reduction of td. This 
indicates that indeed oj^/t^ is the relevant ratio to distinguish the adiabatic from the 
antiadiabatic regime. For oJo/td ^ 1, and independent of oj^/t, the system has physical 
properties that can be described within the simple renormalisation scenario associated 
with the complete Lang-Firsov transformation. 

In figure 11 (b) we make the same observation for A = 3, still with ujq = 1. 
For td > 0.5 the kinetic coefficient L shows the transition familiar to us from the 
previous discussions of intermediate or small phonon frequencies, which is in contrast 
to the physics in the antiadiabatic regime. For <^ 1 a sharp peak occurs in L for 
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Figure 11. (Colour online) All results for t = 1, /i = and = 1- Upper row: 
Kinetic coefficient L as a function of A = £p/2td at A = [panel (a)] and as a 
function of Ep at A = 3 [panel (b)] for different td- In panel (a) L is compared to the 
renormalised non-interacting case (dashed lines, see text). The inset in panel (a) gives 
the corresponding optimal parameter jmin- Lower row: Dot spectral function Addi^) 
and integrated spectral weight S{uj) at A = 3, f,/ = 0.1 for Sp = 2 (c), £p = 3 (d), and 
Sp = A [inset, panel (d)]. 



= A. This is of course the behaviour expected for the antiadiabatic regime, which 
is reminiscent of the non-interacting case for small with a peak of L at A = (cf. 
figure 3). 

The lower panels of figure 11 show the spectral function of the repulsive quantum 
dot at small dot-lead hopping = 0.1 for A > [panel (c)], A ~ [panel (d)], 
and A < [inset panel (d)]. Below the "critical" EP coupling we have 7m,m — 0.2 
and obtain a double-peak structure of A^d because both the first and the third term in 
equation (20) give significant contributions. At Sp = 3, the prefactor of the third term 
vanishes (7mm = 1), and a single-peak structure develops. This polaronic peak is located 
at the Fermi energy and contains all the spectral weight. Therefore L is enlarged more 
than three orders of magnitude. Increasing Sp further the polaronic signal is narrowed 
and shifted away from the Fermi level. As a result L decreases off by five orders of 
magnitude. 

Figure 12 gives more results for the experimentally relevant wide-band case, t > td, 
now in the adiabatic regime. Owing to the values of the parameters /x and ujq, we have a 
situation where ImTiddi^^) 7^ in the whole relevant u region and the spectral function 
is given by equation (28), According to the formula for Im'Eddi^), equation (23), the 
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Figure 12. (Colour online) All results for < = 1, A = 0, id — 0.5, /i = 0, and loq — 0.1. 
Panel (a): £" as a function of 7 for different (inset: minimum of -E as a function of 
Ep). Panel (b): L, jmim and n as functions of Sp. For Sp ~ 0.46 a crossover takes place. 
Panels (c) and (d) Spectral functions Addi^j-!) and integrated spectral weight S{lli) for 
coupling strengths below and above e^. 
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L (a) and particle density on the dot n (b) as functions of Sp for varying chemical 
potential fi. 



cases Ep = 0.1 and Ep < Ep [panel (c)] show the predominance of the first term. The shift 
of spectral weight to negative u becomes apparent for Ep = 0.45, indicating the influence 
of the EP interaction. The spectral functions for Ep > Ep and Ep = 1 [panel (d)] make 
evident the suppression of the first-term contribution and the multi-phonon structure 
according to the second term in equation (23). The maxima of the spectral functions 
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Figure 14. (Colour online) Dot spectral function for A = 0, t = 1. <d = 0.5, = 0.1, 
jjL = —1.9 and varying £p. 



are situated near uo = A. Then again the sudden decrease of L at may be understood 
from equation (52) by the sudden change of ^^^^^(O) at 6^. 

Figure 13, for cuq = 0.1, shows the kinetic coefficient L and the particle density on 
the dot, ra, for /i < 0, whereas the dot spectral function is given in figure 14 for yU = — 1.9 
only. Again, we observe an adiabatic transition in L and n. As /i decreases, the critical 
EP-coupling strength moves to larger values, simply because the effective dot level has 
to be lowered by a larger Sp to roughly match the Fermi level. For Sp < Ep [panel (a)], 
the spectrum lies in the interval [—W — fi,W — fi] ^ [—0.1,3.9], where the first term 
of equation (23) contributes, with the apparent infiuence of EP coupling. For Sp > 
[panel (b)], the spectral weight is shifted to a pronounced peak below the Fermi level. 
Because ujq <W — we find no intervals where haT^ddiyj) = 0. 

Finally, we monitor for the wide-band case the transition induced by an increasing 
dot level A (see figure 15 at small (cuq = 0.1, left-hand column) and intermediate-to-large 
(cuo = 1; right-hand column) phonon frequencies). In both cases Im J^ddi^^) 7^ for all u. 
As discussed above, for A 7^ 0, the maximum in Add occurs near u = A = A— £^7(2— 7). 
In particular, for A = —2 (panel (d)) the spectrum consists practically only of one peak 
at about A — Sp with relatively small linewidth. Hence only a weakly damped localised 
state of the current carrier on the dot exists, having an energy lowering equal to Ep. 
The transition from a localised to a delocalised carrier is accompanied by the shift of 
spectral weight to larger frequencies and the infiuence of the first term in equation (23) 
is recovered. The change of ^^^(0) with A leads to the maximum observed for L in 
panel (b). Because the optimal variational parameter 7mm is a continuous function of 
A with a wide range of values (panel (b)), the effective renormalisation of td and A 
depends on the dot level A itself. In contrast to the result for a complete Lang-Firsov 
transformation with fixed 7 = 1 (cf. figure 5 of (Galperin et al. 2006)), we therefore find 
a shift of the maximum of L{A) by less than Sp and L decreases asymmetrically away 
from this point. However, in accordance with (Galperin et al. 2006), we find no phonon 
side band in L{A). 
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Figure 15. (Colour online) Left column: L, 7mm; and n as functions of A for t = 1, 
id — 0.5, wq = 0.1, Bp = 0.4 and /i = [(panel (a)]. At A"^ w —0.24 a crossover 
takes place. Spectral functions for A below [panel (c)] and above A^ [panel (e)]. Right 
column: X, ^min and n as functions of A for td = 0.5, = 1, £p = 1 and /i = [panel 
(b)]. Panels (d) and (f) give the spectral functions Addi'^) and integrated spectral 
weight S{uj) for various A. 



4. Summary 

In this work, we have presented an approach to transport through a vibrating molecular 
quantum dot, which extends a previously developed description for the many-polaron 
problem. The virtue of this approach lies in an incomplete variational Lang-Firsov 
transformation in which the degree of the transformation is determined self-consistently. 
In this way, our approach can describe polaronic effects on transport away from 
the strong-coupling anti-adiabatic regime. Descriptions based on a full Lang-Firsov 
transformed Hamiltonian are, in contrast, restricted to this limit. 
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With our approach we studied the molecular quantum dot in different regimes, from 
weak to strong coupling and small to large phonon frequency. The dot spectral-functions, 
calculated within a second-order equation of motion approach, allow for a detailed 
analysis of the dynamical properties of the quantum dot in dependence of the model 
parameters. Our results show that the use of an incomplete Lang-Firsov transformation 
is essential to capture the physics for all but very large phonon frequencies: in many 
cases the optimal parameter 7 differs significantly from unity. 

The present study is open for extension in several important directions. On the one 
hand, extension to finite voltage bias is necessary. Since our approach is developed in 
the Green function formalism, this extension, e.g. using Keldysh techniques, is possible 
and will be addressed next. On the other hand, our approach correctly captures the 
physics for a large range of possible parameters, but even with an incomplete Lang- 
Firsov transformation one encounters problems at very small phonon frequency. 

In conclusion, the presented work carries over important concepts and ideas well 
known from polaron physics, especially the crucial modification of the Lang-Firsov 
transformation, to the study of vibrating molecular quantum dots. 
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